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We calculate the scattering lengths of Nambu-Goldstone bosons interacting with D mesons in a covariant 
formulation of chiral perturbation theory, which satisfies heavy-quark spin symmetry and analytical properties of 
loop amplitudes. We compare our results with previous studies performed using heavy meson chiral perturbation 
theory and show that recoil corrections are sizable in most cases. 
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I. INTRODUCTION 

In recent years, studies of charmonium and open-charm systems have witnessed a renaissance. This was mainly led by the 
experimental discoveries of various new particles, either the still largely mysterious X, Y, Z particles or the new open-charm 
states. Many of these states cannot be easily understood in conventional quark models without introducing new degrees of 
freedom in addition to their basic qq structure, notably multi-quark components, i.e. qqqq. 

An interesting resonance in this context is the D* (2317) with a mass of 2317.8 ± 0.6 MeV and a small width of several 
MeV (T < 3.8 MeV) [1]. The Z?* (2317) was first observed by the BaBar collaboration in the inclusive D+ir° invariant mass 
distribution from e + e~ annihilation data at energies near 10.6 GeV [2] and later confirmed by Belle [3] and CLEO H. The 
nature of this state has been extensively discussed in the literature B5I 417II . All studies seem to agree that the coupling of the 
D* (2317) to the nearby DK threshold cannot be ignored. From this perspective, it is particularly interesting to note that the 
D* (2317) can be "dynamically" generated in coupled-channel unitary approaches with interaction kernels provided by either 
chiral perturbation theory (%PT) lfl5l[l6Tl or a SU(4) Lagrangian ifnll . Such approaches have provided many interesting results 
in the past few years, e.g., in explaining the nature of some low-lying hadronic states such as the /o(600) and the A(1405) (for a 
comprehensive list of references, see Ref. fl8ll ). In this dynamical picture of the D* (2317), the interactions of DK and coupled 
channels play a decisive role. A quantity that characterizes the strength of such an interaction at low energies is the scattering 
length. Although it cannot measured directly given the short life time of the D mesons, it can nevertheless be studied on the 
lattice lfl9ll20Tl . The s-wave DK scattering lengths have recently been computed by several authors ll2Tll22ll . In Ref. J2lll . the 
calculation was performed using a covariant formulation of xPT up to next-to-leading-order (NLO) and its unitarized version. 
While in Ref. B22I1 . the calculation was performed using the Heavy-Meson xPT (HM^PT) Il23l427ll in the heavy quark limit up 
to next-to-next-to-leading (NNLO) order, where recoil effects have been neglected. The authors cautioned, however, that since 
the D mesons are not heavy enough, recoil corrections may not be small and have to be studied. 

To our knowledge, recoil corrections have so far not been systematically studied in x?T describing the interactions between 
heavy-light mesons and Nambu-Goldstone bosons. They have been, however, studied quite extensively in the one-baryon sector 
with three flavors u, d, and s. There these corrections were found to be fairly large and play an important role in the studies 
of many physical observables j28l - l3Tll . Because of the large baryon mass that does not vanish in the chiral limit, covariant 
baryon xPT often faces the so-called power-counting-breaking (PCB) problem l32ll . This problem has been traditionally dealt 
with using a dual expansion in terms of both p/A x and 1/Mb, where p is a generic small quantity, Mb the baryon mass and 
A x = An the chiral symmetry breaking scale. This is the celebrated heavy baryon x?T 031 l34ll . Though very successful in 
describing many observables, this approach is not covariant and modifies the analyticity of loop amplitudes. More importantly, 
from a practical point of view it suffers from slow convergence, particularly in the 3-flavor sector. To overcome these problems, 
several other approaches to deal with the PCB problem have been proposed, which among others include the infrared (IR) ll35ll 
and extended-on-mass-shell (EOMS) |36i 13711 renormalization schemes. While the IR scheme was found to introduce artificial 
cuts in certain cases (see, e.g., Ref. IB8I0 . the EOMS approach is fully covariant and conserves the analytical properties of loop 
amplitudes. In a series of applications ll28l - [3Tll . it has been shown that the EOMS approach also improves the convergence 
behavior of SU(3) baryon x?T. 

The main purpose of the present work is to study the scattering lengths of Nambu-Goldstone bosons (0) interacting with D 
mesons in a covariant formulation of xPT by using the EOMS scheme to remove the PCB terms induced by the large D meson 
masses. Given the fact that the D meson mass (~ 1.9 GeV) is much larger than the nucleon mass, it is anticipated that the recoil 
corrections should be smaller than those in the nucleon case. Nevertheless, such recoil corrections may still be sizable as we will 
show in this work. 

This paper is organized as follows. In Section II, we introduce the relevant effective Lagrangians and explain briefly the 
EOMS renormalization scheme. In Section III, we show the numerical results, compare them with those of earlier studies by 
paying special attention to the recoil corrections, followed by a short summary in Section IV. 



II. THEORETICAL FRAMEWORK 

In this section we introduce the chiral Lagrangians relevant to the present study up to NNLO and explain briefly the EOMS 
renormalization scheme used to remove the PCB pieces appearing in the covariant loop calculation. 



A. Chiral Lagrangians 

To introduce the chiral effective Lagrangians, one must specify a power counting rule. In the present case the light meson 
masses and the field gradients d^cj) are counted as 0(p), while d^P, d^P*, rap, and mp. are counted as 0(1), where 
<j> denotes the Nambu-Goldstone bosons, P = (£>°, D+, D+) and P* = (D*°,D*+, D* s + )^ are the D and D* meson fields. 
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The Nambu-Goldstone boson propagator,^3 — T , is counted as Otp 2 ), while the heavy-light pseudoscalar and vector meson 

q 

propagators, g2 _' m ^ and q2 _ l m 2 ^ {-g^ + are counted as ©(p" 1 ). 

The leading order chiral Lagrangian describing the self-interaction of Nambu-Goldstone bosons has the standard form: 

£ (2) = 7 ^(((d t ,m-^) 2 + M<P i ), (1) 
To 



where $ collects the pseudoscalar octet fields 



+ -% tt+ K+ 

V2 V6 



* = ^2| 7T- -^ + ^ A- 1. (2) 

= diag(r7i 2 , m 2 , 2m^ — m\) is the mass matrix, and /o is the pseudoscalar decay constant in the chiral limit. Here and in 
the following, (• • • ) always denotes the trace in the respective flavor space. 

The lowest order chiral Lagrangian for the heavy-light pseudoscalar and vector mesons isQ 



= (V^PV^P^) - m 2 D (PP^) - (D/*"^') + m 2 D (P* u P*^ 

9 

2m 



+i 9 (p; u ^ - p^p; f > + ^({p;u a dpp? - d p p; Ua p^)e^) o) 



where V^Pa = d^P a - T b ^P b and V^P\ = d^P\ + T^P* with a (b) the SU(3) flavor index, g is the heavy-light pseudoscalar- 
vector coupling constant of dimension 1 and m is the mass of the heavy-light meson in the chiral limit. The vector and axial- 
vector currents, and u M , are defined as 

IV = ~(u + d, J ,u + ud IJ ,u + ) and = i(v) d^u - ud^u r ) (4) 



with u 2 = U = exp (^j^j ■ The numerical value of g can be fixed by reproducing the D* + — > D°ir + decay width. Using the 

PDG average, T D , + = 96 ± 22 keV and BR D ,+^ DQ7T + = (67.7 ± 0.5)% [1], one obtains T D ,+^ D07r + = = 

65 ± 15 keV and accordingly g = 1177 ± 137 MeV. 
The NLO Lagrangian relevant to our study is 

C {2) = -2 [c <PP t )<x+) - ^(Px+P 1 ) - c 2 (PP^)(u^) - c 3 (Pu^u^)] , (5) 



where \± = u'Aiu' ± uAiu. Here we have adopted a convention consistent with that of Ref. [22] for the purpose of 
later comparison. In general, there are more terms contributing to 4>P scattering, e.g., (d^Pd^P^) (u^u u ), (d^Pu^u 1 ' '<9„P^) , 
(d^Pu^u^dvP^) . However, at <fiP threshold, it can be easily shown that these terms lead to the same structure as those propor- 
tional to C2 and C3 and therefore can be neglected. The low-energy-constant (LEC) c\ can be determined from the mass splitting 
between strange and non-strange heavy-light mesons within the same doublet, i.e., 

- 8ci(m 2 K - ml) = (m 2 Ds - m 2 D + m 2 D , - mf>„)/2. (6) 

Using the masses given in TableU we obtain c\ = —0.225. 



TABLE I. Numerical values of (isospin-averaged) masses Q|] and decay constants (3^1 (in units of MeV) used in the present study. The 
fx/f-K ratio is consistent with the latest determination |40] while the /,,/ f w ratio is in agreement with that determined in a number of other 
approaches, see, e.g., Ref. |41]. The eta meson mass is calculated using the Gell-Mann-Okubo mass relation: = {^m 2 K — rn^)/3. 



m.D 






m Ds 


m D 




m K 


m r} 


u 


!k 


U 


1972.1 


2112.3 


2008.6 


1968.5 


1867.2 


138.0 


495.6 


566.7 


92.4 


1.22/, 


1.31/ w 



1 Because of heavy-quark spin symmetry, the pseudoscalar and vector mesons can be assigned to the same multiplet. 

2 D* mesons in NLO and NNLO Lagrangians do not contribute to the NNLO D<\> scattering lengths and therefore will not be explicitly shown in this work. 



(») (b) (c) (d) (e) 

FIG. 1. Tree-level contributions at LO [(a),(b),(c)], NLO (d) and NNLO (e). 



At NNLO one has0 

£(3) = -L K (d»P(x-)p1 - P(x-)d»P*) + ^<9"pr M pt - PT ll & i pt){x + ) 

+71 {d»Px+T»P ] - PT^x+d^P^) + 72 {V^PP^ - PV^P^) (2? mX +> , (7) 

where x~ — [x_,it p ]. Although the number of relevant LECs at NNLO is considerably smaller than that present in meson- 
baryon scattering processes ll42ll43ll . it is still relatively large considering the scarcity of lattice data. Therefore in the present 
study we follow the approach adopted in Refs. l44T,l45ll in studies of meson-baryon scatterings and put the NNLO LECs to zero. 
This is acceptable in the present work because we focus on the recoil corrections and not on the absolute values of the scattering 
lengths. 



B. Power counting restoration and the EOMS renormalization scheme 



In a covariant formulation of %PT describing the interactions between heavy-light mesons and Nambu-Goldstone bosons, one 
has to face the PCB problem. That is to say, in the calculation of a loop diagram one may find terms with a chiral order lower 
than that determined by the naive power-counting as prescribed in the previous sub-section. Such analytical PCB terms can 
be removed, just as in baryon xPT, by using the heavy-meson expansion, the IR or the EOMS renormalization prescriptions. 
The essence of the EOMS approach lies in the fact that x?T, by construction, contains all the structures allowed by symmetry. 
Therefore, the PCB pieces appearing in a loop calculation can always be removed by redefining the corresponding LECs. This 
is equivalent to removing the finite PCB pieces directly from the loop results. In practice, this can be achieved in two slightly 
different ways: 1) one can first perform the loop calculation analytically, and then remove the PCB terms, or 2) one can first 
perform an expansion in terms of the inverse heavy-meson mass, l/m#, calculate the PCB terms and then subtract them from 
the full results. It should be noticed that the second prescription is different from the heavy-meson (baryon) expansion because in 
general integration and expansion may not commute. But since the PCB terms are finite and analytical, the second prescription 
should always work. 

In the present study, we have explicitly checked that all the PCB terms appearing in our loop calculation can be removed by 
redefining the LECs introduced in the previous sub-section. 



(a) (b) (c) (d) 



FIG. 2. NNLO loop contributions that survive in the infinite heavy-meson mass (ran — ► oo) limit. 



3 In Ref. l22ll only one such term, that proportional to k, was considered. 
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FIG. 3. NNLO loop contributions that vanish in the mjj — > oo limit. Their crossed counterparts are not shown here but included in the 
calculation. 

III. RESULTS AND DISCUSSION 

The tree-level diagrams contributing to the (dimensionless) threshold T-matrices at LO [Fig. (la,lb,lc)], NLO [Fig. (Id)] and 
NNLO [Fig. (le)] are shown in Fig.Q] Among the three LO diagrams, the contact term [Fig. (la)] yields the following results 

T« = -^, T«=0, T« = 0, = 7f>=0, T« = -™, (8) 

Jo Jo Jo 



_,(!) _ 2m n m D (!) _ (i) _ rh D m K (i) _ m D m K „(i) _ m D m K 

7 ~ f2 ' 1 8 — U ' 1 9 — f 2 ' J 10 — j- 2 ' J H ~~ f2 

J JO ./0 ^0 



where 1, • • • 11 denote the D S K, DK(1), D s ir, DK(0), D s -q, L>tt(3/2), Dir (1/2), Dr], D S K, DK(1) and DK(0) channels, 
respectively. In labeling the 1 1 channels, we have explicitly shown their isospin in the parentheses whenever necessary. In the 
above results, fo is the Nambu-Goldstone boson decay constant in the chiral limit. The contributions given by diagrams (b) and 
(c) of Fig. 1 also count as 0(p), but at threshold they are in fact 0(p 2 ) and have the same structure as those provided by the 
NLO Lagrangians and therefore can be effectively taken into account by redefining the LECs C\ and Co (see below). Hence we 
will not calculate them explicitly in the present work. Following the same arguments, we will also neglect similar diagrams at 
NNLO with one of the LO vertices replaced by a NLO vertex. 

Lagrangian (0 provides the NLO tree-level contributions [Fig. (ld)]0: 



,(2) = Cmi 2 K (2) _ (Co + C^mj, (2 ) _ (Cp + C^ml (2) _ (3Ci~C )to 2 

fo 2 ' 2 Vl ' 3 Vl ' 4 Vl 



K 



T (2) _ (TCi - Co) mg - 16ci (mg - m 2 ) (2) _ drnj (2) _ Cim 2 

°/o Jo Jo 

(2) _ 4ci (m 2 , - m 2 ) + (Cp + 2Ci) to 2 (2) _ fim^ (2) _ G\m| ^ (2) _ Cpm 2 - 

J 8 — of2 ' J 9 — f2 ' J 10 — f2 ' J ll ~~ f2 ' KVL > 

°J0 JO JO JO 

where we have introduced two combinations of the 4 LECs: C\ = 4(2cq — c\ + 2c2 + C3) and Co = 4(2co + c\ + 2c2 — C3). 

The NNLO loop contributions can be separated into two groups: those that survive in the infinite heavy-meson mass (m# — > 
00) limit (Fig. [2]i and those that vanish in the ra^ — > 00 limit (Fig. [3)- For the first group, our results recover those of Ref. ll22ll 
in the tuh — > 00 limit. For the second group, our calculation shows that they indeed vanish in the m# — > 00 limit but are 
not negligible in a covariant calculation, as shown below. It should be noted that all the loop contributions can be calculated 
analytically except the box-diagrams [Fig. (31) and its crossed counterpart]. However, the analytical results are quite involved 



4 These results are the same as those of Ref. [22] except that there the expression for is incorrect 14^1 . 
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and therefore we refrain from showing them explicitly. As explained earlier, in the present work we are going to neglect all the 
NNLO counter terms. Accordingly we have removed from our loop results all the NNLO analytical terms. Therefore our loop 
results contain only NNLO non-analytical and higher-order terms. This way the differences between our covariant loop results 
and those of HMxPT are strictly recoil corrections. 

The s-wave scattering lengths are related to the threshold T-matrix elements through 

a = —, ; rTthr. = ^—7 7 rTthr., (13) 

sir (m^ + nip) aTr(m<j, + mp) 

where is the Nambu-Goldstone mass induced by explicit chiral symmetry breaking and nip is the D meson mass of a 
particular channel. It should be noted that although in the calculation of T t h we have used the average D meson mass, in 
evaluating the scattering lengths [Eq. (fT3l l we use the physical masses as given in Table U We also use the physical decay 
constants instead of the chiral limit ones since the differences between them are of higher order. Furthermore, we have introduced 
T with dimension of a length for a more transparent comparison with the HMxPT results of Ref. 12211 . 

The unknown low-energy constants, C\ and Co, could in principle be fixed by reproducing data, which are however not yet 
available. As in Ref. [22], one can fix them by reproducing the preliminary lattice QCD data (in units of fm): 00-^(3/2) = 
-0.16(4), o D x(i) = -0.23(4), a D ^ = 0.00(1), and a DsK = -0.31(2). 

In Table [TT1 we list the T-matrix elements order by order (T^ for LO and for NLO) and the scattering lengths for the 
1 1 independent (strangeness, isospin) channels computed in the mu — > 00 limit. Since we have neglected the NNLO counter 

~ (3) 

terms, we have denoted the NNLO results by T£ ' stressing the fact that they contain only loop contributions. We have fixed the 
LECs C\ and Co by a least-squares fit of the four lattice points [20]. We obtain a x 2 /dof = 5.0 and the following values for the 
LECs: C\ = 0.4 ± 1.2 and Co = 9.6 ± 10.4 at the 95% confidence level. Clearly, the four lattice data do not constrain well the 
two LECs. 

In TableHni we show the values of the T-matrix elements and scattering lengths found in the covariant approach. The fit yields 
a x 2 /dof — 4.5 and the following values for the two LECs: C\ — 2.0 ± 1.2 and Co = 4.0 ± 10.4. Comparing the relativistic 
and non-relativistic results, one can easily see that the recoil corrections are sizable. For instance, T£ for D S K changes from 
— 1.9 fm to —5.1 fm, and T^ for D s rj changes from 0.1 + 9.7i fm to —4.4 + 5.8* fm when going from HMxPT to covariant 
XPT. 

In obtaining the numbers shown in Tables [TT1 and Hill we have removed the ultraviolet divergences0by the modified minimal 
subtraction (MS) renormalization scheme and have set the renormalization scale fi to 4^/^ . In the covariant x?T, there are two 
heavy scales, A x and nip,. We have checked that our results would remain qualitatively the same if we had set /1 to nip,. 

In Table IIV1 the loop contributions T^ 3 ' calculated in covariant %PT are decomposed into three parts: part A comes from 
diagrams Figs. (2a-2f); part B comes from diagrams Figs. (2g-21); part C comes from diagrams Fig. 3. It is clear that diagrams 
Figs. (2a-2f) provide the most important contributions, while those from diagrams Fig. 3 are similar in size to those from 
diagrams Figs. (2g-21), which is different in HMxPT, where the contributions from diagrams Fig. 3 vanish. 



TABLE II. Threshold T-matrix elements T and scattering lengths a (in units of fm) in the non- relativistic yPT up to NNLO (f (1) for LO, 
T^ 2 ' for NLO and T' 3 ' for NNLO including only loop contributions). The preliminary lattice QCD results [20] are denoted by a(lQCD) and 
have been fitted to fix the two LECs Ci and Co. 



(s,i) 


(2,1/2) 


(1,1) 






(1,0) 


(0,3/2) 




(0,1/2) 




(-1,1) 


(-1,0) 


Channels 


D S K 


DK 


D s -K 


DK 


D sV 


Dn 


Dn 


Di) 


D a K 


DK 


DK 


fm 


-7.7 








15.4 





-3.2 


6.4 





7.7 


-7.7 


7.7 


f(2) 


0.8 


9.7 


1.1 


-8.0 


-1.2 


0.1 


0.1 


7.0 


0.8 


0.8 


18.5 


f.(3) 
1 L 


-1.9 


-1.6 + 5.7i 


-1.1 


3.5 


0.1 + 9.7i 


-0.8 


0.3 


0.8 + 4.8i 


0.2 + 8.5i 


-3.4 


4.8 


a 


-0.28 


0.27 + 0.19i 


0.00 


0.36 


-0.04 + 0.30i 


-0.15 


0.26 


0.25 + 0.16i 0.28 + 0.27i 


-0.34 


1.03 


a(lQCD) 


-0.31(2) 




0.00(1) 






-0.16(4) 








-0.23(4) 





It should be pointed out that since we have neglected all the NNLO counter-term contributions, we are not in a position 
to comment on the convergence behavior of either the covariant or the HMxPT results. Furthermore, because of the nearby 
resonance £>* (2317), a pure x?T calculation, such as ours, in the DK (0) and D s i] channels should be taken with care, where 
coupled-channel unitarity may play an important role (for a relevant discussion, see Ref. J2l|]). 



5 They appear in both HMxPT (see also Ref. |23| 1 and in covariant xPT- The HMxPT results could have been made renormalization scale independent if we 
had kept the NNLO counter terms. 
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TABLE III. Same as Tablelllll but with the threshold T-matrix elements and scattering lengths calculated in relativistic xPT. 



(S.I) 


(2, 1/2) 


(1,1) 






(1,0) 


(0,3/2) 




(0,1/2) 


(-1,1) 


(-1,0) 


Channels 


D S K 


DK 


D S TY 


DK 


D s r) 


Dty 


Dty 


Drj D S K 


DK 


DK 


fW 


-7.7 








15.4 





-3.2 


6.4 


7.7 


-7.7 


7.7 


f(2) 


3.9 


5.8 


0.7 


1.9 


4.9 


0.4 


0.4 


5.2 3.9 


3.9 


7.8 


rp(3) 

1 L 


-5.1 


-2.1 + 4.9i 


-0.7 


2.5 


-4.4 + 5. 8i 


-0.6 


0.3 


-0.4 + 3. 8i -0.9 + 4.4i 


-6.2 


7.7 


a 


-0.28 


0.12 + 0.16i 


0.00 


0.66 


0.02 + 0.18i 


-0.13 


0.28 


0.16 + 0.12i 0.34 + 0.14i 


-0.33 


0.77 


a(lQCD) 


-0.31(2) 




0.00(1) 






-0.16(4) 






-0.23(4) 





TABLE IV. Decomposition of the relativistic NNLO threshold T-matrix elements (P arl; A from diagrams Figs. (2a-2f); part B from 
diagrams Figs. (2g-21); part C from diagrams Fig. 3). 



(s,i) 


(2,1/2) 


(1,1) 






(1,0) 


(0,3/2) 




(0, 1/2) 




(-1,1) 


(-1,0) 


Channels 


D S K 


DK 


D s ir 


DK 


D S 7J 


Dn 


Dtt 


Dr\ 


D 3 K 


DK 


DK 


A 


-6.6 


-2.1 + 4. 2i 


-0.6 


0.7 


-3.5 + 6.7i 


-0.8 


0.5 


-1.8 + 3.4i 


-1.9 + 6.3i 


-7.2 


6.7 


B 


1.9 


-0.6 


-0.1 


1.9 


0.2 


0.2 


-0.2 


0.7 


0.7 


1.5 


-1.8 


C 


-0.4 


0.7 + 0.8i 


0.1 


0.0 


-1.0-0.9i 


0.0 


0.0 


0.7 + 0.5i 


0.3- 1M 


-0.5 


2.8 


f[ 3) = A + B + C 


-5.1 


-2.1 + 4.9i 


-0.7 


2.5 


-4.4 + 5. Si 


-0.6 


0.3 


-0.4 + 3.8i 


-0.9 + 4.4i 


-6.2 


7.7 



IV. SUMMARY 

We have studied the scattering lengths of Nambu-Goldstone bosons interacting with D mesons using a covariant formulation 
of xPT. In particular, we have studied the recoil corrections by comparing the relativistic with the non-relativistic results. Our 
studies show that the recoil corrections are sizable, which should be kept in mind in future studies and in using the HMxPT 
results. Based on available information we cannot conclude which framework is better, although in principle one should trust 
more covariant results, particularly when recoil corrections are large. 

Up to now, xPT describing the interactions between heavy-light mesons and Nambu-Goldstone bosons has often been used 
in the non-relativistic limit. With more precise data and lattice QCD results becoming available, one may have to study more 
carefully the effects of recoil corrections. The present work should be seen as a first step in this direction. 
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